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Abstract 

The article proposes generalizations of the macroscopic model of plasma 
of scalar charged particles to the cases of inter-particle interaction with 
multiple scalar fields and negative effective masses of these particles. The 
model is based on the microscopic dynamics of a particle at presence of 
scalar fields. The theory is managed to be generalized naturally having 
strictly reviewed a series of its key positions depending on the sign of 
particle masses. Thereby, it is possible to remove the artificial restriction 
contradicting the more fundamental principle of action functional addi¬ 
tivity. 

This work was founded by the subsidy allocated to Kazan Federal University 
for the state assignment in the sphere of scientific activities. 

1 Introduction 

In the previous articles the Author considered the statistical systems of scalar 
charged particles and constructed the cosmological models based on such sys¬ 
tems [DlHll]. Particularly, in nig it was obtained the self-consistent set of 
equations describing a statistical ^stem of particles with scalar interaction. In 
the recent Author’s works [5J IHl the macroscopic theory of the statistical sys¬ 
tems has been improved significantly in part of formalism development and also 
it has been extended to the case of fantom scalar fieldfl In [lOl [T^ a question 
about possibility to include negative effective particle masses in the theory has 
been raised. The negative answer to this question has been given in m since 
negative values have been obtained for particle number density in the statistical 
theory at given assumption. However as later more in-depth researches have 
shown, this answer was obtained incorrectly. It turned to be that it is required 
to review a series of key points of the relativistic kinetic theory which depend 
on a sign of particle effective masses in order to strictly resolve the question of 
possibility to include negative effective masses into the theory. These key points 

^see also monographes ElIS]. 

^see also m and review m- 
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appear at earlier stages of the relativistic kinetic theory preceding the defini¬ 
tion of the macroscopic moments. The generalization of the theory to the case 
of multiple scalar fields is being considered in this article. When an attempt 
to make such a generalization is made, it turns out that the developed theory 
contradicts to the fundamental principle of the action function additivity. The 
assumption of non-negativity of particles’ effective mass is proved to be a the¬ 
ory statement which exactly leads to the contradiction. Thus, there raises a 
requirement to revise the kinetic theory of systems of scalar charged particles. 

The results of this revision are presented in the article. 

2 The Dynamics of Particles With The Scalar 
Interaction 

2.1 The Canonical Equations of Motion 

A normalized invariant differential of volume of relativistic particle’s T 8-dimensional 
phase space, being a vector bundle T = X x P with Riemann base X{g) and 
vector bundle P{X), with regard to a pair of canonically conjugated dynamic 
variables a;* (confi^rational coordinates) and Pi (coordinates of the generalized 
momentum) is [13)4 


dr = -^dXdP = 
(27r)3 


dx^ dx^ dx^ dx'^ dPidP2dP3dPi, 

(27r)'^ 

(1) 

where 


dX = ^/^dx^dx^dx^dx*; 

(2) 

dP = _ dPi dP2 dPz dP4 

(3) 

are invariant differentials of volumes of the configurational and momentum 
spaces correspondingly and p is a degeneration factor; for particles with spin S 
it is p = 25'-|-1. Further to shorthand we introduce like phase space coordinates 
rja, a = 1,8: 

r]i = x\ r]i+4 = Pi, (i = l,4), (4) 

where the expression for differential of volume of the phase space ([T]) takes the 
simplest form: 

^ ^ a=l 

(5) 

The canonical equations of relativistic particle motion in 
have the following form (see e.g., [2]): 

the phase space T 


^Here and further the universal system of units is used G = c = h = 1. In ordinary units 
it is 277 —>• 277h. 
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dx'- dH dPi dH 

ds dPi ’ ds dx'^ ’ 

where H{x,P) is a relativistically invariant Hamilton function and u® = dx'‘ jds 
is a particle velocity vector. 

In consequence of antisymmetry of the canonical equations of motion ([6]) and 
symmetry of the phase space o the differential relation [13] known as Liouville 
theorem in classical dynamics, is fulfilled 



(7) 


According to this relation the phase volume of the world tube of particles is 
constant. 

Calculating the full derivative of dynamic variables function 5'(a;®,Pfc) and 
taking into account ([61) we obtain: 


( 8 ) 

where the invariant Poisson brackets are introduced: 

dH dH d'i’ 

Let us note that Poisson bracket (|9]) can be rewritten in the explicitly covariant 
form using the operator of covariant Cartan differentiation, Vi, 0, (see e.g., 

= + ( 10 ) 

where Vi is an operator of covariant Ricci differentiation and P^^ are Christoffel 
symbols of the second kind with respect to metrics gij of base X. Operator V 
is defined in such a way that 

V,Pfc = 0 (11) 

and the following symbolic rule of functions differentiation is fulhlled: 


V,«'(a:,P) = V,[«'(a:],P). (12) 

This rule means that in order to calculate Cartan derivative of function (a:, P) 
it would be enough to calculate its ordinary covariant derivative as if the mo¬ 
mentum vector was covariant constant. Because of this equality the introduced 
operator is quite convenient for execution of differential and integral operations 
in the phase space P. Thus let us write the Poisson bracket (jH]) in the explicitly 
covariant form: 

dH ~ d^ ~ 

Covariant derivative in a bundle P m 

^Cartan covariant derivatives were first introduced into relativistic statistics by A.A. Vlasov 

m- 


3 



Further, as a result of (0 the Hamilton function is an integral of particle motion: 

JTLT 

-j- = [H, H\=Q,^ H = Const. (14) 

The relation ( 11 ^ 1 ) can be called a normalization ratio. Due to linearity of the 
Poisson bracket any continuously differentiable function f{H) is a Hamilton 
function as well. There is a sole possibility to introduce the invariant Hamilton 
function, square by the particle generalized momentum at presence of gravita¬ 
tional and scalar fields. It is the following one: 

H{x,P) = ]^[il)[x){P,P)-ip{x)\, (15) 

where (a, h) here and further is a scalar product of the momentums a and b with 
respect to the base metrics: 


(a, 6) = gika'-b^, 

and 4’{x) and '•p{x) are certain scalar functions of the scalar potentials. Let us 
choose a zero normalization of the Hamilton function in the relation (HI 0112]: 

H{x,P) = ]^[ip{x){P,P) - ip[x)]=Q, (16) 

where from we find: 

ip,p) = r. ( 17 ) 

from the first group of the canonical equations of motion m we obtain the 
relation between the generalized momentum and particle velocity vector: 

u^ = — = 'ipp^^P^ = (18) 

as 

Substituting (I18|) into the normalization ratio dnl), we obtain: 

{u, u) = iptp. 

Hence for particle velocity vector’s normalization ratio fulfillment 

{u, u) = 1. (19) 


it should be: 


Ipip = 1 => "0 = (/?' 


Thus particle’s Hamilton function can be defined by a single scalar function 
^p{x). Taking into account the last relation, let us write the Hamilton function 
in the final form: 


H{x,P) = ^[ip ^{x){P,P)-ip{x)\=Q. 


( 20 ) 
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From the canonical equations (jS]) we can obtain the relation between the gen¬ 
eralized momentum and particle’s velocity vector: 


= ip 


ds 


( 21 ) 


From (I17|) it follows that the vector of generalized momentum if timelike: 

(P, P) = ip^^ 0. (22) 


Let us note a relation that could be useful further, being a consequence of m, 
(fTSll and (HD): 

[P,P'=] = = (23) 

where V® = k is a covariant derivative’s symbol. 


2.2 The Equations of Motion In The Lagrange Definition 

Let us obtain using (ID the equations of motion in the Lagrange defintion [5]: 
d?x'' : dx^ dx^ 


ds^ 


ds ds 


ik 


where: 


(24) 


■p^k ^ pM ^ gik _ 

is a tensor of orthogonal projection on the direction u so that: 


(25) 


(26) 


P^'^Uk = 0; = 3. 

The strict consequence of the velocity and acceleration vectors orthogonality 
follows from these relations and the Lagrange equations 


gtkU —7— = 0 . 

as 


(27) 


Let us notice that the Lagrange equations of motion (1241) are invariant with 
respect to a sign of the scalar function ip{x): 

ip{x) —>■ —ip{x). (28) 

The Hamilton function dini) at its zero normalization is also invariant to the 
transformation (1^51) . Therefore from (HD, HU) and the Lagrange-Filer equations 
(1^ it follows that tp scalar’s square has a meaning of a square of the effective 
inert mass of a particle, m* in the scalar field: 


2 2 
ip =m^. 


(29) 

Let us notice that the following action function corresponds to the cited choice 
of the Hamilton function: 


S= m^ds. 


(30) 


It formally coincides with the Lagrange function of the relativistic particle with 
a rest mass m* in the gravitational field (see e.g., [H])- 
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2.3 The Choice of A Mass Function 


Let there be now n various scalar fields <l>r in the system and each particle has 
n corresponding fundamental scalar charges qr {r = l,n) among which there 
can be some zero charges. Then the problem of choice of a function 
appears. Not specifying this function for now, let us note the next important 
circumstance: For equations of motion ([6]) to allow the linear integral of motion 
= (^, P) = Const, according to ([8|) it is enough and sufficient that [H, 'L] = 0. 
This, in turn, is possible if and only if [lOj : 

(^, P) = Const L ipgik = 0, (31) 

where L is a Lie derivative in the direction ^ (see e.g., [T5]L 

Let us consider statistical fields gik and <i>r allowing timelike Killing vector 
= SI when the total energy of charged particle is conserved P 4 = Eq = 
Const > 0. Let us choose a system of coordinates where ga 4 = 0,a,f3 = 1,3 
so that coincides with the world time t. Then from the relations connecting 
the kinematic velocity vector u® with the total momentum of a particle Pi m 
it follows: 

P4ds = tpdt. 

Therefore, if we require to conserve the same orientation of the world time t and 
proper time s (i.e. = dt/ds > 0) then the mass function should be chosen in 

such a way that it is always nonnegative: 

m* = \q}\ > 0. 

However such a conservative approach had been being used in Author’s papers 
DP - 0 and seeming at first sight the correct one, turns to be contradicting 
to the more fundamental principle of the Lagrange function’s additivity. As 
it was shown in [10], the negativeness of the particle’s effective mass function 
does not lead to any contradictions at the level of microscopic dynamics since 
the observable momentum of particle (as well as the 3-dimensional velocity 
v°‘ = u°^/u^) conserves its orientation as opposed to the unobservable kinematic 
4-velocity of a particle u*: 

p-=rn,^ = P\ (32) 

ds 

Then the choice of a linear mass function in pOI) corresponds to the principle 
of action function’s additivity: 

TO* = mo + ^ (33) 

r 

where toq is a certain initial rest mass and is a charge of particle relative 
to the scalar field $*. which we presume functional independent. At toq = 0 the 
entire rest mass of particles is ensured by interaction with scalar fields: 

m^, = (fi = ''^qr^r- (34) 
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This choice satisfies also aesthetic criteria since in this case the Hamilton func- 
tion lfTSl) does not depend on the rest mass. From the other hand it is clear that if 
choose function </?('&) in form (|55)) the Lagrange equations (1241) become symmet¬ 
ric with respect to the change —>■ —or qr —>■ —qr- Therefore if particles 
a and antiparticles d are different in sign of scalar charges then m* = — m*, 
but trajectories of particle and antiparticle in gravitational and scalar fields are 
identical. In such case according to (ITSl) 4-vectors of their kinematic velocity 
are different in sign: u® = — u®, and have coinciding vectors of the observable 
generalized momentum P® = P®. Let us note that as this takes place, vectors of 
the 3-dimensional velocity {a = 1, 3) for particles and antiparticles 

also coincide with each other: = u“. Nevertheless we do not assume mo = 0 

for now postponing this for the future. 

Having chosen the effective mass in form (l33l) . the Hamilton function (|20l) 
and the normalization ratio (1221) for the generalized momentum take form: 


Hix, P) = ^ [m* ^ix){P, P) - m*] = 0, 


(35) 


(P,P)=m^ (36) 

Let us notice identities valid for the Hamilton function (1551) which could be 

useful in future: 

VjP = -Vim*, (37) 

[P,4-] = —P®V,4' + a,m*|^, (38) 

m* aPi 

where 4'(a;,P) is an arbitrary function. 


2.4 The Quantum Equations 

From the classical Hamilton function (1551) with a help of standard substitution 

I 

Pi —>■ ih\7i (39) 

we obtain the Hamilton operator: 

H = -^m*-i(?iV"V,Vfe-|-m^). (40) 

Thus for a free massive scalar field we obtain wave equations in form of standard 
Klein-Gordon equations with the only difference being bosons rest mass changed 
to the effective mass m- 

(□-f m*/ft^)4 = 0. (41) 

For free fermions we obtain corresponding Dirac equations: 

(S7®V,-fm*)4 = 0, (42) 

®In this section we temporarily renounce the use of the universal system of units where 
h= 1 . 
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where 7 are spinors. 

Let us notice that from (1411) at substitution 'L = 4> and choice of the mass 
function m* = as well as for the mass function m = |( 7 <I>| it straight away 
follows the equation of the scalar field with cubic nonlinearity: 

□4> + = 0. (43) 


Thus the self-action constant in the equation of the scalar field takes well-defined 
meaning: 



i.e., it is defined identically to the fine structure constant for the electromagnetic 
field. 


3 Calculation of the macroscopic means 

In order to extend the general relativistic kinetic theory to the case of negative 
effective masses of elementary particles it is required to review the key points 
of the relativistic kinetic theory which could depend on a sign of particle mass. 
It has been revealed that there are only two of such key points and they are 
defined by two circumstances: the relation between proper microscopic time of 
a particle and proper time of macroscopic observers and the relation between 
particle velocity vector and its generalized momentum. 

3.1 The Invariant Distribution Function 

Generalized formalism of invariant distribution functions was developed in arti¬ 
cles [SITS]. It is required to thoroughly apply this formalism to the case being 
considered in order to account the possibility of the negative sign of the effective 
mass of particles. According to this formalism, to define the macroscopic means 
in the relativistic phase space it is required to define the unit timelike field of the 
macroscopic observers Ui{x) : (U, U) = 1. The synchronization of measurement 
acts for individual particles is carried out by the clocks of such observers. This 
timelike field in its turn defines certain spacelike three-diemensional surface V 3 , 
displacements along which (5a;* are orthogonal to this field: 

F 3 : Sx^Ui = 0, (44) 

and displacements along this field dx"^ define the synchronized proper time r of 
observers: 

^ = c/* = l^dT = dx^U,. (45) 

dr dr 

Thus in the observers’ macroscopic reference frame: 


X = VxT^dX = dVdr. 


( 46 ) 


Apparently, the connection between the proper time of particle s and the 
synchronized proper time r of observers in each point of the configurational 
space is established by the relation: 


dr dx^ 
ds * ds 


(47) 


Now taking into account relation (15^ being the consequence of the canonical 
equations (0), we finally obtain: 


dr 

ds 


= ^{U, P) => 


1 ds 
m^:{s) dr 


1 

{U{T),Pis))- 


(48) 


The relation (1491) can be considered as a differential equation in the function 
s(t), resolving which we define the relation between proper time of particle and 
time measured by the clocks of synchronized observers: 


s = s(r). 


(49) 


The invariant 8-dimensional identical particles distribution function F{x, P) 
is introduced in the following way Let phase trajectory of particle being 

defined by the canonical equations m in the phase space T = P{X) x X is: 

x" = x'-{s)-, Pj = Pj(s) ^ ? 7 a = ? 7 a(s) (a = 1,8), (50) 


where s is a proper time of particle. Then the number of particles registered by 
the observers in range dT of the phase space can be defined as [81II3]: 

dN{T) = F{x, P)S{s - s{T))dr. (51) 


Let us notice that particle number is a scalar depending, however, on the choisce 
of observers’ field t/*, i.e., on the choice of reference frame in Riemann space 
X, while 8-dimensional distribution function itself F{x,P) being introduced 
by the relation m, is an invariant in phase space T. Let us note also that 
it is impossible to provide any other definition for the invariant distribution 
function in the 8-dimensional phase space. All definitions of this function which 
were introduced earlier are special cases (ISTl) realized in the dedicated reference 
frames. 


3.2 The Macroscopic Means of Dynamic Functions 

The definition of the invariant distribution function (ED and other dynamic 
functions together with it is a first key point of the relativistic kinetic theory 
which depends on the sign of particles’ mass. Let 'tlj{x,P) = il’i'd) is a certain 
scalar function of dynamic variables. Then according to (IHTll its macroscopic 
mean 4'(r) in range C T is defined in the following way: 

4'(t) = J 'tlj{r]{s))dN = 

n 

J Fir]{s))ilj{r]{s))S{s - s{T))dT. (52) 

n 
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Let us emphasize one more time that distribution function F(x, P) is an invari¬ 
ant in the phase space as opposed to the macroscopic means which are always 
defined with respect to observers’ field U. Putting further in accordance with 
(Hill , (HSl) X = V xT ^ dX = dV dt, let us write expression (15^ in the explicit 
form: 

V T 

f dP'ip{r))F{T])'il^{r))S{s - s(t)) (53) 

P{X) 


To carry out integration over t in (j53l) let us take into account the relation 
between t{s) (15^ and r(s) (HHl) and the property of Dirac d-function: 


6{s — s(r))dt = 


dr 

ds 


S(t — T)ds 


= TO* 


-1 


{U, P)S(t — T)dt. 


(54) 


At derivation of (iMl) we took into account the fact that generalized momentum’s 
orientation does not depend on the sign of the effective mass as opposed to 
kinematic velocity vector’s orientation. Taking now the account of ((5^ and 
((551) and conducting integration over time coordinate, we finally obtain: 


4'(t) 


2S+1 f 

(27r)3 J 


UdV 


1 

| to *| 


X 


j P^dPtlj{r])F{r])ip{ri) 

P(X) 


(55) 


Particularly, supposing '0 = 1, we obtain the total number of particles in range 
P: 


NiV) 


25'-bW 

7^ J 


UdV 


1 

| to *| 


X 


J P^Fir])dP = J {U,n)dV, 

P{X) V 


(56) 


where particle flux density vector is introduced: 


o c I 1 r 

n\x) = 37 / P^F{r^)dP 

(2'7r)3|TO*| J 

P{X) 


(57) 


This factor of effective mass modulus was not taken into account in the previous 
articles because of above cited conservatism reason. 
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3.3 Transformation to the 7-Dimensional Distribution Func¬ 
tion 

As a result of generalized momentum normalization ratio (1361) the invariant 8 - 
dimensional distribution function F{x, P) is a singular one on the mass surface. 
Therefore it is necessary to introduce a distribution function non-singular on 
this surface f{x,P) by means of the relation: 

F{x,P)= 6 {H{x,P))f{x,P). (58) 


Calculating the relation 


Fix, P)dP ^-=dPidP 2 dP 3 dP 45 (H(x, P))f(x, P) 

with a help of Dirac 5-function’s properties and the explicit form of the Hamilton 
function (I55|) . we find: 

F{x, P)dP = 

\m 




dP4dP2dP3'-^5{P4 - P:)f{x, P) X 


n 


= |TO*|(iPo^(C4 - P4)f{x,P), 


where P^ = P 4 is a positive root of the normalization equation (I36|) and 


dPa = 


1 dPidP2dP3 

7^ 


(59) 


(60) 


is a differential of volume of the 3-dimensional momentum space. 

As a result formulas for the macroscopic means (I55p can be written through 
the 7-dimensional distribution function f{x,P) in the following form: 

^(r) = J U.dV J P^dPoijiv)fiv), (61) 

V P+(X) 

where it is necessary to substitute P 4 by the positive root of the mass surface 
equation and P+ is an upper part of the mass surface pseudosphere. Thus at 
transformation to the 7-dimensional distribution function the explicit depen¬ 
dency on the effective mass disappears. 

Therefore the next symbolic rule being understood in a sense of integration 
over corresponding phase volumes, is valid: 

'0(?7)P(?7)5(s — s{T))dr —>• 

Hv)fmU,P)dVdPo, (62) 

where 77 are dynamic variables on the 6 -dimensional subspace T 0 (r) = V x Pq G 
T with the differential of volume: 


dTo = dVdPo. 


(63) 
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In particular, for the particle number flux density ?;ecto)0 (1641) from (ICTl) we find: 

2S' + 1 f 

= J P^f{ri)dP. (64) 

Po(X) 

Let us notice that in m the incorrect conclusion was made regarding the 
absence of dependency between particle number flux density vector and the 
effective mass sign. This conclusion was a result of the fact that macroscopic 
means’ calculation has been started with a calculation of distribution function 
moments while in case of negative masses the correct calculations should have 
been started on the earlier stage, namely they should have been started with 
integral relations of form (1521) . 
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